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DISCOVERY OF AN EXTENDED DEBRIS DISK AROUND THE F2V STAR HD 15745 
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ABSTRACT 

Using the Advanced Camera for Surveys aboard the Hubble Space Telescope, we have discovered 
dust-scattered light from the debris disk surrounding the F2V star HD 15745. The circumstellar disk 
is detected between 2.0" and 7.5" radius, corresponding to 128—480 AU radius. The circumstellar disk 
morphology is asymmetric about the star, resembling a fan, and consistent with forward scattering 
grains in an optically thin disk with an inclination of ^67° to our line of sight. The spectral energy 
distribution and scattered light morphology can be approximated with a model disk composed of sili- 
cate grains between 60 and 450 AU radius, with a total dust mass of 10 -7 M© (0.03 M$) representing 
a narrow grain size distribution (1—10 /im). Galactic space motions are similar to the Castor Moving 
Group with an age of ~10 8 yr, although future work is required to determine the age of HD 15745 
using other indicators. 

Subject headings: stars: individual(HD 15745) - circumstellar matter 
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1. INTRODUCTION 

The collisional debris from solid bodies orbiting main 
sequence stars results in circumstellar dust disks that 
produce excess emission at thermal infrared wavelengths 
(|Backman fc Parescd Il993). Excess thermal emission 
may also arise from a star heating interstellar dust 
(e.g.. iKalas et al.l I2002T ) and therefore obtaining images 
of dust-scattered light is a key method for constrain- 
ing the origin of dust. Furthermore, direct images of 
dust-scattered light may be used to obtain the color 
and polarization properties of the grains, which depend 
on t heir size distribution , composition and structure 
(e.g. iGraham et alj 120071 ). The overall geometry of a 
debris disk may also reveal sources of dynamical per- 
turbations, such as those from ste l lar and sub-stellar 
companions (e.g.lRoaues et alJll994HLiou fe Zooklll999t 
iMoro-Martin fc Malhotrall2002j ). 

Here we show the first scattered-light images of a de- 
bris disk surround ing HD 15745, an F2V star at 64 pc. 
ISilverstond (|2000f) first reported HD 15745 as a debris 
disk candidate star based on excess thermal emission 
Li R /L+ — 1.8 x 10~ 3 using IRAS data, and Ltr/L^ = 
1.3 x 10~ 3 using ISO data. IZuckerman fc Sons) (|2004l ) 
modeled the spectral energy distribution (SED) with 
a single temperature dust belt at 27 AU radius with 
L IR /L* = 1.2 x 10~ 3 . 

2. OBSERVATIONS & DATA ANALYSIS 

We observed HD 15745 on 2006 July 17 with the Hub- 
ble Space Telescope (HST) Advanced Camera for Sur- 
veys (ACS) High Resolution Channel (HRC) corona- 
graph. The star was acquired behind the 1.8" diameter 
occulting spot and three 700 second integrations with 
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the F606W (A c = 591 nm, AA = 234 nm) broadband fil- 
ter were obtained at a fixed position angle. The three 
frames were median combined for cosmic ray rejection. 
The point spread function (PSF) was then subtracted us- 
ing images from the HST archive of five F stars obtained 
with the ACS coronagraph at the same wavelength and 
comparable integration times. All five versions of a PSF- 
subtracted image of HD 15745 show the fan-shaped cir- 
cumstellar nebulosity. The five versions were then me- 
dian combined to produce a final image that was then 
corrected for geometric distortion. We also detected HD 
15745 from Gemini North with the MICHELLE cam- 
era on 2006 September 20 using the N' filter. We mea- 
sure a flux of N'= 74.7±2.9 mJy within a 1.4" radius 
aperture, with an aperture correction based on near- 
contemporaneous images of the photometric standard 
star (HD 20893). 

3. RESULTS 

The PSF-subtracted HST data reveal a fan-shaped 
nebulosity surrounding HD 15745 in the approximate re- 
gion between P.A.=190° and P.A.=10° (Fig. 1). We 
interpret the nebulosity as an azimuthally symmetric cir- 
cumstellar disk inclined to our line of sight and composed 
of dust that has an asymmetric scatter ing phase func- 
tion t hat produces a fan-like appearance ([Kalas fc Jewittl 
1996). Another possibility — asymmetric grain erosion 
and blowback by encounters wi th the interstellar medium 
(|Artvmowicz fc Cla mpin 1997])— is less likely given that 
the southeastward proper motion (^t a =46 mas/yr, n$=- 
47 mas/yr) and locatio n of the star well within the lo- 
cal interstellar bubble (|Kalas et al.ll2002h . although at 
least one exce ption has been found for the latter point 
(jGaspar et al.ll2007f ). 

In principle the detection limit for the inner radius of 
the nebulosity is set by the occulting spot radius (0.9"), 
but the presence of quasi-static speckles produce signifi- 
cant PSF subtraction artifacts, limiting the inner detec- 
tion radius to ^2" (128 AU) radius. The outer radius 
of the disk is a sensitivity-limited value. In data that 
have been binned 8x8 pixels (0.2" x 0.2"), the disk is 
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detected as far as ~7.5" (480 AU) radius. The 23.0 mag 
arcsec -2 surface brightness isophote (Fig. 2) shows a rel- 
atively symmetric fan-like morphology. The radial sur- 
face brightness profiles (Fig. 3) are monotonic and can be 
approximated by power laws with an index of —3.3 ±0.1. 

4. DISK MODELS 

We model the SED and scattered light appearance 
of the HD 15745 disk using a Monte Carlo radiation 
transfer code (MCFOST; Pinte et al. 2006) in which 
a star radiates isotropically in space and illuminates an 
azimuthally-symmetric parameterized disk. In a first 
step, wc determine the stellar parameters by fitting the 
SE P up to 11 Mm (Fi g. 4). Using NextGen mod- 
els (|Baraffe et al.l Il998) for the stellar photosphere, we 
find that the best fit model has T e // = 6800 K and 
R+ = 1.31-Rq. 

We then proceed to fit the entire SED of HD 15745 
to determine the disk properties. In order to con- 
sider as few free parameters as possible, we adopt a 
p(r, z) ex r a exp(— z 2 /2<r 2 ) number density distribution, 
with a = — 1 and a constant <r z = 2AU, that extends 
from R in to R out = 450 AU and sums up to a total disk 
mass Md- The choice of a z — 2 AU is consistent with 
disk heights in ferred from the edge-on debris disks Pi c 
and AU Mic dKalas fe Jewittl Il995l iKrist et al.l l200l . 
and does not significantly influence our results for values 
o~ z > 1 AU (which is required to keep the disk optically 
thin to starlight). The r^ 1 number density profile fits 
the observed slope of the surface brightness (see below). 
We also fix the dust grain properties to be a collection 
of com pact spherical silicate grains from iDraine fc Led 
(1984), with size distribution N(a) oc a -3 5 from a m i„ to 
a>max — 10a mjn . Because the disk is assumed to be gas- 
free, each grain size has its own temperature: smaller 
grains are hotter than larger at a given distance from 
the star. Therefore, for each value of a m i n , we need to 
adjust Ri n to obtain a maximum dust temperature in 
the range 60-80 K, i.e., a peak dust emission around 40- 
50 /im (Fig. 4). Once this is set, we can adjust the total 
dust mass to fit the SED. For values of a m i n ranging from 
0.125 to 8 /im (with a factor of 2 increment), we perform 
this 2-parameter adjustment using a grid of 9 regularly- 
sampled values of Ri n and Md and compute the associ- 
ated x 2 value. The narrowness of the disk thermal emis- 
sion peak requires a max < 10 /im (hence a m j n < 1/im). 
Larger grains deviate from the Rayleigh regime in the 
mid-infrared and result in a much broader emission peak 
than is observed. By allowing a smaller value of a m i n 
and larger values of a ma x/ 'o-mini we do not find a lower 
limit to a m in from the analysis of the SED alone, down 
to at least 0.1 /im. Provided Ri n and Md are adjusted 
accordingly, we always find satisfying fits to the SED. 

We then explore the disk surface brightness in order to 
better constrain a m i n . Indeed, assuming a max /a rn i n = 10, 
the scattering properties of the dust populations con- 
sidered here vary strongly with a m i n - The albedo at 
0.606 /im decreases from 0.82 to 0.65 as a min increases 
from 0.125 /nm to 1 /im, because of the depletion of the 
highly reflective small grains. At the same time, the scat- 
tering asymmetry parameter, g — (cos#), increases from 
0.65 to 0.89 because large grains scatter preferentially 
forward. As a consequence, the models with small values 
of a m in can be rejected based on our ACS image: sim- 



ulated images of the disk with a m i n = 0.125 /im are ^4 
mag arcsec~ 2 too bright compared to our observations. 
Furthermore, with such a dust population, the predicted 
flux ratio between the front and back side is low enough, 
independently of the inclination, that it should have been 
detected in our data set. To fit both the surface bright- 
ness of the disk and the apparent front/back asymmetry, 
our model indicates that a m i n > 0.5 /im, implying a nar- 
row grain size distribution and justifying a posteriori our 
assumption of a max /a mm = 10 

Figure|4] illustrates our best fitting model, which has 
a m in = 1/im, a max = 10 /im, R in = 60 AU, M d = 
10~ 7 M Q and L/a/L*=2.2x10~ 3 . Assuming a disk in- 
clination i = 65 — 70°, a power law fit to the modeled 
surface brightness profiles results in a ~-3.5 index and 
the observed front/back intensity ratio at 2'.'3 along the 
disk semi-minor axis is > 10. Both figures are consistent 
with the observations. Subsequent models were calcu- 
lated with a m i n — 0.75 and 1.25 /im to better constrain 
this parameter. We obtained marginally acceptable sur- 
face brightness fits in both cases and we conclude that 
0.75 /im < a m i n < 1.25 /im in the disk surrounding HD 
15745. This result is consistent with a grain blowout size 
of ~1 /im calculated for silicate grains. Associated uncer- 
tainties on the disk inner radius and total disk mass are 
on the order of 10%. Our quantitative estimates of a m i n 
and a max depend on the assumptions we made, but the 
overall conclusion (order of magnitudes for grain sizes, 
hence inner radius and total dust mass, and the need for 
a narrow size distribution) are robust. 

5. DISCUSSION 

HD 15745 is now the fifth main sequence F star 
with a debris disk resolved in scattered light. Three 
of these have ring-like structure, whereas HD 15745 
shares a wide-dis k architecture with HD 15115 (F2V; 
iKalas et al.ll2007i ). The three ri ng-like architectur es are 
found around H D 139664 (F5V; IKalas eTaIl l200^. HD 
18132 7 (F5/F6: iSchneider et alJ 12006ft. and HD 10647 
(F9V; IStapelfeldt et alJl2007l) ~At the present time, the 
dust properties have been modeled only for HD 181327, 
with a Henyey-Gr eenstein scattering asy mmetry parame- 
ter g = 0.3±0.03 ^Schneider et al.ll2006f ). For debris disks 
around stars of other spectral types, AU Mic (M1V) has 
a scatt ering asymmetry par ameter similar to that of HD 
15745 ([Graham et al.ll2007t ). wh ereas most other d ebris 
disk, such as Fo malhaut (A3V: IKalas et al.l 120051 ) and 
HD 107146 (G5V; lArdila et al.ll2004h have alow scatter- 
ing asymmetry, comparable to that of HD 181327. Our 
present model suggests a concentration of dust grains in 
the range 1 — 10 /im, which is consistent with a wavy grain 
size distribution (departure from a simple power law; 
e.g. Campo Bagatin et al.l Il994t iThebault &: Augereaul 
2007]). Small grains do not exist in sufficient quantity to 
counteract the highly asymmetric scattering phase func- 
tion of the larger grains, whereas in other debris disk sys- 
tems the collisional cascade may be robust to the smallest 
sizes, producing relatively symmetric scattering. Further 
refinements for the properties of dust around HD 15745 
and other debris disks requires resolved imaging obser- 
vations at other wavelengths and with polarization. Ul- 
timately the explanation for why one debris disk has a 
particular grain size distribution may depend on differ- 
ences in disk dynamical properties and the composition 
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of grains (e.g., volatile vs. refractory elements). 

Among the F stars discussed above, HD 15115 and 
HD 181327 may bel ong to the Beta Pic M oving Group, 
with age -12 Myr (|Zuckerman et al.ll200lh . HD 15745, 
on the other hand, may be a member of the Cas- 
tor moving group bas ed on its Galactic space mo- 
tion. iMoor et alJ 1200(1 calculate (U, V, W) = (-16.5 ± 
1.1, -10.8±1.3, -10.7±0.7) km s" 1 for HD 15745, consis- 
tent with {U, V, W) = (-10.7±3.5, -8.0±2.4, -9.7±0.7) 
km s" 1 given for the Ca stor Moving Group (CMG; 
iBarrado v Navascuesl [T998) . Among the CMG mem- 
bers, HD 38678 (( Lep), another well-known debris 
disk star has a Galactic space velocity (U, V, W) = 
(-14.4±3.6, -11.1 ±3.0, - 8.4± 1.8) km s" 1 that i s clos - 
est to that of HD 15745. IBarrado v Navascuesl (|1998h 
estimate that the CMG has an ag e of 200±100 Myr, 
although IZuckerman fc Song] (|2004l ) estimate a younger 
age for HD 15745 ("30?" Myr) based on its location 
on an A-star H— R diagram and the relatively large dust 
optical depth. Moor et al. (2006) similarly conclude that 
nearly all stars with Lm/L* > 5 x 10~ 4 are younger than 
10 s yr. The infrared excess for HD 15745 is an order 
of magnitude greater than that of CMG members Fo- 
malhaut, Vega, and £ Lep. From the five CMG can- 
didates that are F-stars, two have have been detected 
at far-infrared wavelengths. HD 119124 (F7.7V) has 
Lir/L* = 3— 6xl0~ 5 , compar able to that of Fomalhaut 
and Vega dTrilling et al.1 12001. H P 130819 (F3V) has 
L m /L* = 3 x 10 4 (jChen et ai~ll2005[ L which is still a 
factor of seven less than that of HD 15745. 

Recent work demonstrat es a large spread in dust op- 
tical depths at a given age ^R hee et a 1 2007) which may 
be a consequence of the spread of initial disk masses 
(|Wvatt et al.ll2007h . Therefore the high optical depth of 
HD 15745 cannot uniquely discriminate between a ~10 
Myr or a ~100 Myr age. Assuming steady-state evo- 
lution, the dust optical depth of HD 15745 could arise 
from a 10 Myr-old debris disk that started with ^10 
M ffi of solid material, or a 100 Myr-old disk that started 
with ~100 Mq of solid material. The ambiguity is aug- 
mented by allowing stochastic spikes in dust production 
that may occur in a few per cent of debris disks, (e.g. 
iBeichman et al.ll2006| ). A more precise age for HD 15745 
therefore requires further input using stellar properties. 
However, we note that the large geometric vertical op- 
tical depth of our model (r = 5.8 x 10~ 3 at 60 AU ra- 
dius) implies a fiducial collisional destruction timescale 
t c ~ 1/ Sir ~ 10 5 yr at 60 AU. This is significantly shorter 



than the age estimates for HD 15745, implying that the 
circumstellar dust is collisionally evolved. Moreover, the 
rapid collision timescale indicates that the radius of the 
inner disk hole approximates the location of the planetes- 
imal belt that regenerates the observed dust, instead of 
signifying the radius where grains inspiraling from Poynt- 
ing Robertso n drag are trapped in resonances with plan- 
etary bodies (jWvatt et al ' 
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6. SUMMARY 

Using the H ST ACS coronagraph, we present the first 
images of dust-scattered light from a debris disk sur- 
rounding the F2V star HD 15745. The fan-shaped mor- 
phology is consistent with an azimuthally symmetric disk 
with asymmetric isophotes due to an asymmetric dust 
scattering phase function combined with a 67° inclina- 
tion to the line of sight. HD 15745 belongs in a category 
of wide disks with a sensitivity- limited outer extent near 
500 AU. Our silicate grain models for the SED and the 
disk surface brightness suggest an inner disk radius of 60 
AU, a somewhat narrow grain size distribution between 1 
and 10 /xm, and a total disk mass of 10 -7 M Q . Future ob- 
servations should determine the optical and near-infrared 
colors, as well as polarization, to further constrain grain 
properties. High-resolution imaging of the inner part of 
the disk should reveal the inner hole that in the model 
presented here has radius 60 AU, or 0.5") radius from the 
star projected on the sky. Due to the asymmetric scat- 
tering phase function, the inner boundary of the disk 
hole should appear as a 1" x 2" elliptical feature. Our 
silicate-grain disk model predicts a peak U-band and H- 
band surface brightnesses of 14.8 and 12.7 mag arcsec -2 , 
respectively, or about 1 magnitude redder than the star 
itself. Near-infrared data would thereby test our present 
disk model that is dominated by relatively large grains 
(>1 /im), where the albedo is not significantly dimin- 
ished in the near-infrared, yet the near-infrared scatter- 
ing is more isotropic (g=0.57 in H vs. 0.89 in V) as the 
observing wavelength approaches the grain size. 
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Fig. 1. — False-color, PSF-subtracted image of HD 15745. North is up, east is left. Observations were made using the F606W filter 
and the Advanced Camera for Surveys High Resolution Channel. Gray fields cover the focal-plane masks in the ACS HRC. The axis of 
symmetry for the nebulosity is PA~280° . 




Fig. 2. — Surface brightness isocontours for the HD 15745 debris disk converted from F606W to Johnson V-band (derived using 
STSDAS/SYNPHOT with a Kurucz model atmosphere and the appropriate observatory parameters). North is up and east is left. The 
outermost contour is 23.0 mag arcsec -2 , with a contour interval of 0.5 mag arcsec -2 . 
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Fig. 3. — V-band radial surface brightness cuts between 2.5" and 7.0" radius along the northern (PA=360°) and southern (PA=210°) 
portions of the HD 15745 nebulosity. We extracted photometry in boxes with sizes of 0.1" (along the radial direction) by 0.5", and we 
plot the median value from the five subtractions of the HD 15745 PSF. Representative error bars show the standard deviation of these five 
values for each box. The noise due to PSF subtraction artifacts is greater for the south side of the disk, as evidenced by the rippled form 
of the radial profile on scales of a few tenths of an arcsecond and the larger error bars. In some radial regions the south side of the disk 
may be fainter than the north side by ~0.4 mag arcsec -2 , but these surface brightness asymmetries should be verified in future data that 
have a higher signal-to-noise ratio. 
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Fig. 4. — Observed SED of HD 15745 along with our best-fitting model (see text for model parameters). Our new N' (A c =11.2 /im, 
AA=2.4/xm) Gemini measurement is indicated with a diamond (Gemini program ID: GN-2006B-C-12). HD 15745 is unresolved with a 
full-width at half-maximum of 0.3" in the N-S direction, and 0.5" in the E-W direction. Photometric calibration was accomplished by 
observing HD 20893 in N' and assumi ng an extinction of 0.172 magnitude per air mass. The 24 — 90 fim data represent the IRAS and 
ISO values given by Moor et alj pOOlih . The dotted and dashed curves represent the stellar and disk emission, respectively. The inset is a 
log-stretch image of our best-fit model, observed with an inclination of 67° and rotated so as to match the orientation of the disk in Fig. 
1. The white circle represents the 2" radius region that is hidden in Fig. 1 and 2. 



